Arbitrary-order lensless ghost imaging with thermal Hght 
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Arbitrary A/'th-order {N > 2) lensless ghost imaging with thermal light has been performed for the 
first time by only recording the intensities in two optical paths. It is shown that the image visibility 
can be dramatically enhanced as the order N increases. It is also found that longer integration 
times are required for higher-order correlation measurements as N increases, due to the increased 
fluctuations of higher-order intensity correlation functions. 
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Compared with the first "ghost" imaging experiment 
with two-photon entangled light second-order ghost 
imaging and ghost interference with thermal light |2|, 0, 
i, la i, 0, i, U EO, [HI has a low visibility which theo- 
retically can never exceed 1/3, and in fact will be much 
lower than 1/3 in practical applications. Moreover, when 
a 2-dimensional high resolution image of a complex ob- 
ject is required, the better the resolution, the worse will 
its visibility be. This is one of the limitations for the 
practical application of thermal ghost imagin g (G I) . For- 
tunately, recent studies [H, [H, [H, [H, |l3, 111 on the 
higher-order intensity correlation effects of thermal light 
show that the visibility can be significantly improved by 
increasing the order N. In this way, the drawback of low 
visibility in correlated imaging with thermal light can be 
overcome. 

Third-order GI with thermal light has been theoreti- 
cally analyzed to a certain extent [H, [l^ , but recently 
Liu et al pointed out that it is inappropriate to assume 
that second-order correlations play the entire or dom- 
inant role [11. In their investigations of higher-order 
thermal ghost imaging and interference Liu et al showed 
that it is A^-photon bunching that characterizes the A/'th- 
order correlation and leads to the high- visibility in A/'th- 
order schemes. The necessary condition for achieving a 
ghost image or interference pattern in A/'th-order inten- 
sity correlation measurements is the synchronous detec- 
tion of the same light field by different reference detec- 
tors. Multi-photon interference experiments have been 
carried out by Agafonov et al [14], verifying the con- 
clusion that the visibility limits of three-photon and 
four-photon interference are respectively 82% and 94% 
for classical coherent light, as predicted theoretically by 
Richter p. 5]. Cao et al discussed A'-th order intensity 
correlation in double-slit ghost interference with thermal 
light and proposed a scheme to study the visibility and 
resolution of the fringes with two detectors [H. How- 
ever, in their actual experiment only one CCD detector 
was employed, and the measurements were taken first 
with and then without the double-slit in place. Similar 



high-order schemes to obtain higher visibility but for GI 
were also suggested by Agafonov et al a little earlier [11 . 

In this paper we report the first demonstration of an 
arbitrarily high A/'th-order lensless GI experiment with 
pseudothermal radiation. We do not actually need A^ 
light paths but measure the A/'th-order intensity corre- 
lations by means of just two detectors. Moreover, when 
certain conditions are met, no lens is required for obtain- 
ing a well- focused image of the object. 

Generally, in an A/'th-order intensity correlation mea- 
surement, the light beam needs to be divided into A^ 
parts, each of which passes through an optical system 
and then is registered by a detector. In the high- 
intensity limit, the normalized A/'th-order intensity cor- 
relation function is given by [13, 
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where Ij {xj ) is the instantaneous intensity at position Xj 
in the transverse direction, and (...) stands for ensem- 
ble averaging, in practice achieved through time averag- 
ing. It is well known that {x) = (/^(x)) / {I{x))^ = 
N\ for polarized thermal light when all the points in 
Eq. [T] are the same space-time point [U, [22|. This 
A'-th order intensity correlation function can thus 
be obtained by measuring the autocorrelation function 
of /(x), as has been done in the 2nd order GI and 
interference-diffraction experiments with pseudothermal 
light [1, fl6l . 111, with mere knowledge of the inten- 
sity distribution I{x) of that space-time point. It is 
thus reasonable to envisage a setup in which we have n 
beams of one intensity distribution (say passing through 
an object) and N — n beams of another distribution 
(as reference arms), so that the instantaneous inten- 
sities Ii{xi) = I2[x2) = ... = In{xn) = I{yi) and 

In+l{Xn+l) ^n+2(^n+2) ... lN-n{xN-n) HV^)- 

Then the transverse normalized A/'th-order correlation 
function may be defined as 
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where the intensity correlation is composed of an n-fold 
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FIG. 1: (Color online) Experimental setup for arbitrary-order 
lensless GI with a pseudot hernial light source. G: rotating 
ground glass plate; BS: beamsplitter; O: object. 



intensity product at position yi and an (A^ — n)-fold prod- 
uct at position ?/2, so there are — 1 ways to measure 
the A/'th-order correlation function [l^, • This equa- 
tion implies that, in GI, an arbitrary order correlation 
function can be obtained with mere knowledge of the 
two intensity distributions I{yi) and /(^2). 

Here we consider a simple experimental scheme such as 
that shown in Fig. 1, which is basically the same as that 
for lensless 2nd-order correlation ghost imaging [ll[. To 
pass to higher-order intensity correlation measurements, 
the signals of the test and reference detectors must be 
raised to the powers n and N — n. It should be men- 
tioned that the correlation function measured this way 
is not the same as that in an A^-port HBT interferome- 
ter, but the difference vanishes in the high-intensity limit 
as classical intensity fluctuations (excess noise) become 
stronger than quantum fluctuations. We consider the in- 
tensity distributions I{yi) and I{y2) of the speckle fields 
on the Zi and Z2 planes located, respectively, at dis- 
tances zi and Z2 from the thermal source. For simplicity 
but without loss of generality, in the derivation below 
the image is considered to be one-dimensional in the Zi 
plane. 

If Tijj) indicates the field transmission function of the 
object, and a bucket detector is used in the test (object) 
arm, the bucket signal S is 



/ 



S= / J(t/i)|T(yi)rdyi. 
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If we assume uniform illumination so that {I{yi)) = /, 
the average value of S is 
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in which Aobj is equal to the object area when |T(^i)| is 
a step function equal to 1 inside the object and outside. 
In the classical limit we can also write the A^th-order 
correlation function for GI as 
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where we have assumed that all the n test arms have the 
same configuration, i.e., contain the same objects in the 



same positions and identical bucket detectors. It is evi- 
dent that this requirement can be simply satisfied with 
a single test beam and the detector output divided into 
n parts. Since the same reasoning can be applied to the 
reference arms, we can thus perform any A^th-order cor- 
relation measurement with a simple experimental setup 
similar to that for 2nd-order lensless imaging in which 
only two light paths are needed, as shown in Fig. [TJ The 
reference detector output is correspondingly divided into 
N — n parts [2^ . From the viewpoint of photodetection, 
the relation between the photocurrent i{x^t) and light in- 
tensity I{x,t) is [m i{x,t) oc (1/T) j^!^^^I{x,t^T)dT, 
where T is the exposure time of the detector. If s and i 
are the photocurrents of CCDl and CCD2, respectively, 
then from Eq. [5] the A/'th-order correlation function for 
GI is 
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This is the expression that is actually used in our data 
processing. The total intensity of the light source is as- 
sumed to be equal to the sum of the intensities of all the 
test and reference paths. 

An outline of the experimental set-up is shown in 
Fig. [TJ A CW laser beam with a wavelength A of 532 nm 
and beam diameter D of 3 mm is projected onto a 
ground-glass plate G rotating at a speed of 0.7 rad/s 
to form pseudothermal light. The scattered beam is sep- 
arated by a 50% — 50% non-polarizing beamsplitter (BS) 
into two beams. The transmitted beam passes through 
the object, and the reflected beam is the reference beam. 
The two beams are detected by the charged-coupled- 
device (CCD) cameras CCDl (Imaging Source DMK 
21BU04) and CCD2 (Imaging Source DMK 31BU03), re- 
spectively. The object (a mask), which is the Chinese 
character for "light" shown in Fig. [21(a), is placed in the 
test arm where CCDl plays the role of a bucket detec- 
tor. The distance zi between the source and the object is 
equal to the distance between the source Z2 and CCD2, 
namely, zi =2^2 =240 mm. Thus the coherence area A at 
the plane Zi is about 16 /im^, according to the well- 
known relationship of A ^ (Xzi/D)'^; zs is the distance 
between the mask and CCDl and is equal to 70 mm for 
the high-order imaging experiments. The cameras are 
operated in the trigger mode and are synchronized by 
the same trigger pulse. The data are acquired with an 
exposure time (0.1 ms) much shorter than the correlation 
time of the light source, which is on the order of 0.2 s, 
and saved through a USB cable to a computer. 

Figures 2(b) and (c) show the images of the mask ob- 
tained by direct exposure on the camera CCDl alone 
for Zs = 20 mm and 70 mm (not in the bucket de- 
tection mode), after averaging over 20,000 frames. We 
see that Fig. 2(c) is a completely blurred image due to 
interference-diffraction effects, while Fig. 2(b) is quite a 
clear projected image since the CCD camera is close be- 
hind the mask. 
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FIG. 2: First row: (a) Object mask, (b) and (c): Projection 
images obtained by CCDl alone, averaged over 20,000 frames, 
for (b) Z2> — 20 mm and (c) zs — 70 mm. Second row: High- 
order ghost images obtained with thermal light in a lensless 
setup, averaged over 140,000 frames, corresponding to (d) 2nd 
order (AT = 2,n = 1), (e) 10th order {N = 10,n = 9), (f) 20th 
order \n = 20, n = 19). 



Though we cannot obtain a first-order image in the 
case of ^3 =70 mm, we can obtain clearly discernible 
higher-order ghost images (TV > 2), as shown in the 
lower half of Fig. 2. For these shots a total of 140,000 
frames was grabbed by each camera per plot at a rate 
of about 1.35 Hz. The arbitrary-order correlation images 
are calculated according to Eq. [6] by a Matlab program 
in which the total intensity s of each frame recorded by 
the bucket detector CCDl is divided by n to give the 
intensity s/n of each component test arm, and similarly 
the output i{y2) of CCD2 which is captured at the same 
instant is divided by TV — n to give the intensity of each 
component reference beam; then the product of (s/n)^ 
and {i{y2)/{N — n))^~^ is averaged over all the exposure 
shots and divided by (s/n)^ {i{y2)/{N — n))^~^ to give 
the normalized A/'th-order correlation function ^^\y2)- 
Here {s) and {i{y2)) are obtained by averaging all the 
frames of CCDl and CCD2, respectively. In our experi- 
ment n = — 1 is chosen, as this gives the best visibil- 
ity with the minimum number of frames to be processed. 
The reconstructed 2nd, 10th and 20th order ghost images 

which are normalized following ^n^\y2)/{ln^\y2))Max 
are shown in Figs. 2(d), (e) and (f), respectively. It can 
be clearly seen that the image visibility increases with 
increasing as predicted. This is one of the chief ad- 
vantages of GI systems. Another advantage is that the 
position distribution of the light intensity transmitted 
through the mask is not necessary in any-order GI, so 
long as all of the light is captured by the bucket detec- 
tor. This has been demonstrated by Meyers et al who 
succeeded in obtaining 2nd-order ghost images even with 
a turbid distortive medium placed between the object 
and bucket detector |24| . 

It should be remarked that the rate of increase of the 
visibility in object imaging is much lower than that in in- 
terference experiments, such as the A^th-order Hanbury 




FIG. 3: 4th-order lensless ghost images for n 
(from left to right). 



1,2, and 3 



Brown-Twiss and interference experiments reported by 
Cao et al [17]. This is because the number Mobj = 
Aobj/Acoh of coherent areas falling inside the object, 
which is defined in Ref. [13], is very high and the visi- 
bility decreases with the increase of Mobj- Moreover, it 
should be noted that there is, in fact, a different number 
of coherent areas M^^^ within the object for each differ- 
ent order intensity correlation measurement if we define 
the coherent area (the resolvable area under the focused 
condition) as the square of the full-width-half-max 
of the A^th-order intensity correlation function. For a 
given A^, M^^^- is also different when n is different. It 
has been shown that the visibility increases with the or- 
der N when N is not very high, while it is independent of 
the resolution and close to unity at very high A^ [13, [3 • 
It should be pointed out that in A/'th-order lensless 
GI there are at least A^" — 1 experimental configurations, 
since n may be different. We take the 4th-order as an 
example to show the difference between them. In Fig. 3 
the 4th-order ghost images for n = 1,2, and 3 are ob- 
tained by processing the same number of frames (i.e. us- 
ing the same integration time) after the exposures with 
the CCD cameras. It is evident that the image is com- 
pletely blurred in the case of n = 1 while a good image 
is obtained when n = 3. It is found that the higher the 
order A^, the longer is the integration time required for a 
clear image to be obtained. Moreover, for a given N the 
lower the n, the longer is the integration time required. 
This can be inferred from Eq. [6] where we can see from 
the self-correlation intensity product i^~^{y2) between 
the reference signals (or for the test signals) that as 
n decreases {or N — n increases), the fluctuations of the 
higher-order correlation functions will increase. Since the 
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FIG. 4: (Golor online) Experimental setup for arbitrary-order 
lensless GI with a pseudothermal light source. G: rotating 
ground glass plate; BS: beamsplitter; M: mirror. 



4 






• g2 




A g4 




• 96 




♦ g8 




* gio 



15 20 25 30 35 40 45 50 55 
No. of pixels 



FIG. 5: (Color online) Experimental results of 2nd, 4tli, 6tli, 
Sth and lOth order lensless GI along one cross-section. Trans- 
verse coordinates are the number of pixels, and each point on 
the plot is an average of the raw data from 5 pixels. The solid 
lines are Gaussian fits. 



fluctuation of the total intensity recorded by the bucket 
detector (the A^th-order self-correlation function of s) is 
less than the sum of the fluctuations of each pixel in the 
reference detector (^(^2)), as we have indeed observed 
experimentally, it is best to choose a large n to obtain 
clearer images and save integration time. 

We have also performed an experiment for a smaller 
Mobj using the experimental setup shown in Fig. [H A 
He-Cd laser of wavelength 441.6 nm is reflected by a mir- 
ror (M) and used with a rotating ground glass plate (G) 
to form the pseudothermal light source, which is about 
1 mm in diameter. A pinhole is placed 18 cm away 
from the ground-glass plate to limit the transverse size 
of the beam, which is separated by a 50% — 50% non- 
polarizing beamsplitter (BS) into two beams. The re- 
flected test beam is totally internally reflected by prisml, 
passes through BS and the object, and is detected on the 
lefthand side of a CCD camera (Imaging Source DMK 
31BU03). The transmitted reference beam is reflected by 
prism2 and then by BS to be recorded on the righthand 
side of the detector. Both prisms can be translated to 



adjust the beam path distances between the light source 
and the detector. The object is a double-slit mask with 
slits of width 150 /im and separation 570 /im. As before, 
for lensless imaging the distances from the source to the 
object and reference detector are equal, =2:2 =354 mm. 
To make Mobj small enough, the data is only collected in 
one dimension along the transverse direction. The exper- 
imental results are shown in Fig. [5l in which n = 
and the 2nd, 4th, 6th, Sth and 10th order image cross- 
sections are obtained by averaging about 4,000 exposure 
frames. We can see that the image visibility increases 
much faster than in our previous experiment, mainly be- 
cause here the number of features Mobj is much smaller. 

In conclusion, we have reported the first demonstration 
of arbitrary-order {N > 2) lensless GI with pseudother- 
mal light by only recording the intensities in two optical 
paths. The experimental results demonstrate that the 
image visibility can be dramatically enhanced as the or- 
der increases. This overcomes the bottleneck of low 
visibility in 2nd-order ghost imaging with thermal radi- 
ation. It is conceivable that higher order GI could also 
be realized with a single detector setup using a laser and 
computer generated spatial modulation [25, 26], but for 
a real thermal light source the field distribution in the 
reference arm must be actually measured by at least one 
nonvirtual detector. As we have shown, in fact just one 
reference detector is sufficient. It is also found that longer 
integration times are needed as TV increases, due to the 
increased fluctuations of higher-order intensity correla- 
tion functions. Moreover, for a given N the integration 
time increases with the decrease of n, so it is best to 
choose a large n to save time. The experimental setup 
for high-order GI is the same as that for the 2nd-order ex- 
periment, and we only require a program to calculate the 
image to any desired order of correlation. This makes GI 
with a thermal source even more promising than before 
for practical applications. 
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